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ABSTRACT 


The feasibility of pulse torquing a pendulous accelerometer 
in a ternary mode is investigated. A comparison is made between 
basic ternary and binary systems. It is shown that a ternary sys- 
tem is basically more stable than a binary system in that stable 
oscillations are smaller. Thus errors in velocity indication are 
less for the ternary system. 

The describing equations of a pulse torqued pendulous ac- 
celerometer are transformed into a non-dimensional form. By so 
doing, general conclusions are reached. Phase plane analysis is 
used to investigate the dynamic performance of the device. An 
error analysis based upon the phase plane analysis follows. 

After conditions necessary for larger limit cycles are inves- 
tigated, a method is devised for selecting system parameters such 
that larger limit cycles are prevented. 

Finally, methods of compensation and error correction using 
Switching logic are investigated. Representative switching logics 


are presented. 
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CHAPTER 1 


INTRODUCTION 


fete Lis TORY 


A basic requirement of any inertial navigation system is that the 
system be able to obtain velocity. One common device developed to 
obtain velocity inertially, is the pendulous accelerometer. In sucha 
device, a torque generator current required to keep a pendulum at 
null, is proportional to acceleration along a sensitive axis. When such 
current is integrated with respect to time, an analog of velocity 1s ob- 
tained. 

With the addition of a digital computer to the inertial navigation 
system, a requirement is generated for an analog to digital converter. 
Such converters are inherently limited in accuracy of conversion, 

As an alternative to providing a converter, a binary pulse torqued 
integrating accelerometer has been under study and development. Such 
an accelerometer has an output which consists of a train of pulses. 
These pulses may be fed directly to a digital computer. Each pulse 
represents an incremental velocity change along the sensitive axis of the 
accelerometer. By summing the train of pulses, taking into considera- 
tion their algebraic sign, a velocity indication is obtained. 

In the operation of a binary pulse torqued integrating accelerom- 
eter, certain errors are introduced which are caused by inherent lag 


in the pendulum servo system. 
1.2. A PENDULOUS ACCELEROMETER 


In general, a pendulous accelerometer consists of a pendulous 
mass suspended in a viscous fluid (Figure 1-1). The pendulous mass is 
constrained to rotate about a single axis. An error detector is used to 
measure the deviation of the pendulous mass from a reference axis. 
After amplification, the error signal is used to command a torque gen- 
erator in such a way that the pendulum is driven back toward its null 
position. The applied torque may, or may not be, a linear function of 


null angle. 
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The equation of motion of the pendulum is 


Pu 
dA dA _ _ 
es aPcosA+ap,, P sin A (1.1) 


J 
where 
P= Pendulosity of the pendutind 
A = Angular position of the pendulum 
J = Moment of inertia of the pendulum 
C = Coefficient of viscous friction 
M,. = Torque Generator torque 


a = Component of acceleration along the input axis of the 


accelerometer 
aR A = Component of acceleration along the Reference Axis 
of the accelerometer 
i = tonnes: 


For the purpose of this report, accelération is defined as the 
second derivative of position with respect to time as viewed in a free- 
fall non-rotating coordinate system. This is not a restriction in that 
Equation (1. 2) is valid in any coordinate system if the acceleration, a, 
is defined to include gravitational attraction and apparent accelerations 
associated with rotating frames. 

The null angle, A, 1s maintained sufficiently small to justify small 


angle approximations being made. With this in mind, Equation (lL. l)~ 


becomes: 
? . 
ee ea -aPt+ap, PA (1. 2) 
2 
dt dt 
Rearranging and integrating: 
t t t 
v- I Ale { my, at-2 94 Say! ap, Adt (1. 3) 
O P Yo ‘8 P dt P O 


where v is the velocity increase from time zero to time t. dA/dt and 
A are assumed zero at time equal to zero, 


By definition: 


-12- 





t 
Vv. = A ( M,_ dt = indicated velocity (1, 4) 
i =p ay) tg 
O 
J ees rea A = null error (1. 5) 
‘pdt 8=6P 
t 
ae \ ap A A dt = cross-coupling error (1. 6) 
e=e +e_ = total error. Clee) 
n Cc 


Substituting (1.5) and (1.6) into (1. 3): 
1 ( 
Vv =— Vi ac ce (1. 8) 
Pp Yo tg n eC 


1.3. A, BINARY PULSED INTEGRATING ACCELEROMETER 


Figure 1-2 is a block diagram of a binary torqued integrating ac- 
celerometer. The pendulum is always torqued by a current of constant 
magnitude. The sense or direction of the current through the torque 
generator may be changed, however. The error detector is sampled 
at arate synchronous with a computer clock. At each sample ~instant, 
the current to the torque generator is switched or not switched as 
necessary to drive the pendulum toward null. Lag is introduced into 
the restoring system by the use of a discrete time between samples. 

The state of the applied torque is based upon the error signal at 
a Sample instant. Practical error detectors have a finite insensitive 
zone near their null position. If at a sample instant, the pendulum is in 
the insensitive zone (hereafter called the dead zone), the applied torque 
is maintained in a sense based upon the last available information. Thus 
additional lag is introduced. 

As a result of lag in the error channel, the pendulum oscillates in 
stable modes or limit cycles. These limit cycles if large allow the null 


error, e€., to be several velocity increments in magnitude. 
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Figure 1-2. A Binary Pulsed Pendulous 
Integrating Accelerometer. 





The chapters to follow investigate the feasibility of reducing the 
errors ina pulse torqued integrating accelerometer by torquing the pen- 


dulum in a ternary mode. 
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CHAPTER 2 


A TERNARY PULSED INTEGRATING ACCELEROMETER 


2.1. GENERAL 


The size of the limit cycle of a pulsed accelerometer may be re- 
duced by one of two methods. First, the lag which causes the ex- 
cessive limit cycles may be reduced. Second, lead may be introduced 
to compensate for lag already in the system. 

If, in a pulse torqued accelerometer, the pendulum null angle is 
such that the acceleration reaction torque is in a direction to reduce the 
null angle, it can be anticipated that the null angle will eventually be 
driven to zero without the use.of additional torque. Thus, the use of 
applied torque to aid acceleration reaction, coupled with the discrete 
time increments at which it must be applied, introduces unnecessary lag 
to the system. This elementary observation prompted this investigation. 

If torquing to aid acceleration reaction torque is to be prevented, 

a ''no torque'' mode of operation must be introduced. Since both posi- 
tive and negative accelerations are to be expected, positive and negative 
torque modes must be retained. Thus a ternary mode of operation is 


suggested. 
2.2. COMPUTER CONSIDERATIONS 


The computer most usually used in inertial guidance systems is 
the Digital Differential Analyzer (DDA). The Digital Differential Analyzer 
can be mechanized to handle binary or ternary inputs. The ternary 
scheme is more accurate in that it has less round-off errors, but re- 
quires more equipment to realize. 

The DDA must accept information ata clock rate. Thus the re- 
quirement that the torque pulses be synchronous with a computer clock 


will be maintained in a ternary system. 
2.3. INFORMATION AVAILABLE 


The accuracy of a pendulous accelerometer depends upon the devia- 
tion of the pendulum from its null position being maintained small. With 


such small deviations allowed, pendulum velocity is extremely havd to 


a= 





measure through system noise. For this reason, it is considered that 
pendulum velocity is not directly available. 
Pendulum position is available if the pendulum is not in the error 
detector dead zone. If the pendulum is in the dead zone, the fact that 
it is in the dead zone is known. This fact may be used in switching logic. 
The past history of applied torque is available and may be useful 
in that it represents the past motions of the pendulum. Since input ac- 
celeration is a slow changing function with respect to pendulum motion, 
the past history of pendulum motion might be useful in predicting future 


motion. Thus, lead might be introduced. 
2.4. POWER SUPPLY CONSIDERATIONS 


If torque pulses are to accurately indicate the same velocity quan- 
tum with repeatibility, the torque generator current must be maintained 
extremely constant. ‘The current taken from the power supply is con- 
stant ina binary system. A ternary system requires that the current 
be switched on and off. Because of the variation in current taken from 
the power supply, a better power supply would be required for the ternary 
system, 

An alternative to switching the power supply off might be to switch 
it to adummy load. At any rate a ternary system would require more 


or better power supply equipment. 
2.5. A BASIC TERNARY PULSED INTEGRATING ACCELEROMETER 


Figure 2-1 shows a basic ternary pulsed integrating accelerometer. 
Switching logic is based upon pendulum position. If at a sample instant 
(clock pulse), the pendulum is not in the error detector dead zone, the 
pendulum is torqued for the next sample interval, ina direction toward 
null. If the pendulum is in the dead zone at a sample instant, the pen- 
dulum is not torqued for the next sample period. 


For a pulsed system, Equation (1.8) can be written: 


(2.1) 


where 
te = time between samples 


7 
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Figure 2-1, A Ternary Pulsed Pendulous Integrating Accelerometer. 
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+1 when torquing in a positive direction 
6. = 40 when not torquing 


-l when torquing in a negative direction 
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The digital error is described by Figure 2-2, It is caused by an 
instantaneous addition of velocity quantum in the computer. The actual 
torque-time pulse has taken place over a finite time interval. 


The total error is now defined by: 


(2, 2) 


D 
UI 
ay 


e 
nt Sa 


2.6. NON-DIMENSIONAL FORM 


It is desirable to express the equations of motion as free of as 
many system constants as possible. By so doing, general conclusions 
can be easier reached. By proper substitutions of relations in Table 2-] 
into Equations (1.2), (2.1), and (2, 2): 


g+o=(6-rtrpy ¢) (2, 3) 
n 
v=) 6é,- E,- E,- Ey (2. 4) 
i=0 
EE +E, +g (2, 5) 
where oo coe 
diz 


Two system‘parameters must be specified before an analysis of 
system performance can be made. Equation (2.3) will have to be piece- 
wise integrated between samples. Therefore 2 ee the ratio of sample 


time to pendulum time constant, must be specified. The non-dimensional 


avon 





Digital Error 
<q 


Figure 2-2, 


Example of Digital Error. 
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TABLE 2-] 


DIMENSIONAL TO NON- DIMENSIONAL TRANSFORMATIONS 


Parameter Dimensional Non-Dimensional Relation 
Notation Notation 
Angle A d + cian 
MJ 
Time t = seek 
J 
Acceleration a r r -aP 
Velocity Vv V eed 
MJ 
Welocity Error e E Id a 
MJ 


BOA 





deac zone, @,, must be specified in order to tell where switching will 


take place, 


2.7. PHASE PLANE ANALYSIS 


In general, the acceleration is essentially constant when viewed 
over several sample intervals. Thus, acceleration is considered con- 
stant in this analysis. For the purpose of analysis, r is assumed much 
greater than TRA @ in prescribing the motion of the pendulum. This does 
not imply that the cross-coupling error is negligible but rather that the 
dynamic behavior of the pendulum is not greatly affected by the cross- 
coupling acceleration torque. This type analysis cannot be made for 
zero acceleration along the sensitive axis in that Tha is the only forcing 
function when the pendulum is not in a torque mode. 


With ther » pa @ approximation, Equation (2.3) becomes: 


¢+¢=(6-1) (2.6) 


Since 6 remains constant between torque switches, Equation (2, 6) 
may be integrated as a first order linear equation in ¢, with results being 


valid between torque switches: 


b=$¢,e 7 +(8-r)l-e7) (2.7) 
where b, is the pendulum velocity at 7 = 0. 
=o +6 (l-e 7) +(8-rl rte 7-1) (2.8) 


where %o is the null angle at 7 = 0. 


It is believed that a phase plane (4 versus @ plot) study lends the 
best understanding of second order nonlinear equations. It is assumed 
that the reader has a basic understanding of phase plane analysis. Since 
an elastic restraint term is absent from Equation (2.6), a phase plane 
template can be constructed for a given value of (6 - r). Such a template 
can be used over and over again by sliding it along the ¢ axis. This is 


possible because the phase trajectory is independent of the relative 
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position of the @ = 0 axis. Time must be plotted as a parameter along 
the trajectory in order to tell when switching will take place. 

Since the null error has been defined as the sum of angle and angle 
rate, it can be found for any point on the phase plane by projecting the 
point in question onto an axis which is perpendicular to all 
b + @ = constant lines. 

Adding Equations (2.7) and (2. 8): 


E = 6 +6, +(6-r)7. (2. 9) 


Equation (2.9) states that the null error is a linear function of time 


between torque switches. 


2.8. PREVENTING APPLIED TORQUE FROM AIDING 

ACCELERATION REACTION TORQUE 

It is desirable that torque not be applied in a direction that will aid 
the acceleration reaction torque. Appendix B shows that if the non- 
dimensional dead zone, PA» is greater than the non-dimensional sample 
ime , TS: a limit cycle cannot exist for zero input acceleration. In gen- 
eral, if PA is greater than ps the pendulum in steady state, will re- 
main out of the switch zone that will switch torque to aid acceleration 
reaction (Figures 2-3 through 2-12). For accelerations greater than 
zero, the forcing function of Equation (2.6) tends to bend the phase tra- 


jectory away from the undesirable zone, 


2.9. EXAMPLE 


In any analysis of a pulsed accelerometer, the non-dimensional 
sample time and the non-dimensional dead zone are functions of system 
parameters. For an illustrative example of analysis by phase plane 
methods, the non-dimensional sample time and the non-dimensional dead 
zone are both taken as unity. The non-dimensional dead zone is assumed 
greater than the non-dimensional sample time by an infinitesimal amount, 
such that limit cycles cannot exist for zero acceleration. 

Figures 2-3 through 2-l2 are phase plane plots of the system re- 
sponse for various values of non-dimensional acceleration. For magni- 


tudes of non-dimensional acceleration equal to or greater than one, the 
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restoring torque cannot overcome the acceleration reaction torque, The 
system becomes unconditionally unstable under such conditions. For ac- 
celerations within the designed range of the accelerometer, the non- 
dimensional acceleration will never have a magnitude greater than one. 
Phase plane plots of a binary torqued system are shown in Figure 
2-3 through 2-12 for comparison purposes. The binary and ternary sys- 
tems have identical system parameters. As could be expected, the limit 
cycles of the ternary system are smaller than those of the binary system 
for the same input accelerations. Higher order limit cycles exist for 
both systems. Such limit cycles give rise to errors in velocity indication, 
For the associated accelerations, all possible limit cycles have been shown 
in Figures 2-3 through 2-12. The conditions under which higher order 
limit cycles may exist in a ternary system are investigated in Chapter 3 


and Appendix C., 


2,10. ERROR ANALYSIS 


From phase plane studies, it is noted that limit cycles for a given 
input acceleration have a certain freedom in position about the switching 
line, The exact amount of freedom that a simple limit cycle possesses 
has been determined and plotted as a function of Te in Appendix C. 

Figures 2-13 through 2-22 are plots of null error versus time for 
the phase plane plots of Figures 2-3 through 2-12, Again, a plot is shown 
for a binary system for comparison. The error plots are not unique in 
that the initial conditions determine the exact location of the limit cycle. 
The null error is bounded by a region determined by the freedom in po- 
Sition of limit cycle, Pes plotted in Appendix C. 

It is noted that the peak nuil error is always less for the ternary 
system, usually by as much as one to two velocity counts. 

The cross-coupling error is equal to the product of the component 
of acceleration along the reference axis, and the average value of the null 
angle. Both systems are biased in that the average null angle is not zero 
for accelerations greater than zero. In general, the average null angle 
is greater for the ternary system if the acceleration along the input axis 
is small and is greater for the binary system if the acceleration along 


the input axis is large, The bias in null angle is equivalent to the forced 
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dynamic error ina pendulous accelerometer with linear feedback. 

The maximum magnitude of the average null angle is always less 
than the magnitude of the dead zone. The maximum value of the cross- 
coupling error is 


E 


(2,10) 








=o d 
ae RA "A 


In general, the digital error, E,, is of opposite sign than the 


null error, thereby reducing the otal oror. The digital error is 
the error introduced by an incremental summation being made in the 
counter, instead of a true integration of torque. 

The error equations used in this analysis have assumed that the 
pendulum is at rest in its null position when the problem commences, 
An exact analysis will require that the error associated with initial 
conditions be subtracted from the final error to determine the actual 


error in a particular problem. 


2,11. REDUCTION IN ERRORS 


It is apparent that errors could be reduced if the size of the limit 
cycle is reduced. Two approaches will be taken to this problem. First, 
system parameters will be adjusted in an attempt to minimize the limit 
cycle. Second, a system of compensation by switching logic will be 


investigated. 
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CHAPTER 3 


MINIMIZING LIMIT CYCLE SIZE BY PROPER SELECTION 
OF SYSTEM PARAMETERS 


3.1, GENERAL 


In the previous chapter it was shown that a ternary pulsed pen- 
dulous accelerometer may operate in several modes of oscillation about 
a bias position. The bias, or average pendulum position, produces two 
distinct errors. First, the bias itself represents an error in velocity 
indication in that a certain amount of acceleration reaction is Opposed 
by viscous friction in the process of the pendulum moving from its null 
position. This error is non-accumulative and its effect may even be 
negligible. Second, cross-coupling error is introduced by the bias 
angle. For small angle approximations, the cross-coupling error is 
the time integral of the product of the bias angle and the component of 
acceleration along the reference axis. Cross-coupling error is ac- 
cumulative (i.e. it may grow with time). It may be minimized by main- 
taining the bias angle small, or may be "computed out" if the bias angle 
is known. 

The additional error associated with the oscillation of the pendu- 
lum about its bias position, is oscillatory and therefore non-accumulative. 
The instantaneous error caused by this oscillation might even be negli- 
gible. However, these oscillations if large may cause the velocity in- 
formation to be presented to the computer in an undesirable form. 
Additional errors may be generated by the computer which are ac- 
cumulative. The nature and extent of these errors will not be investi- 
gated here. To minimize these errors it is desirable to maintain the 
pendulum in its lowest mode of oscillation. (A limit cycle ‘3 of the 
lowest mode isdefined to be that limit cycle which encloses the least 
area of the non-dimensional phase plane. Thus for a non-dimensional 
acceleration of 2/7, a1:2:1:3 limit cycle is of lower mode than a 2:5 or 
me l:1:4 limit cycle. ) 

It is presupposed that certain characteristics are specified before 


initial design of a ternary pulsed accelerometer begins. Sample time, 
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maximum acceleration and physical size and weight limitations may be 
among the specifications. The designer is at liberty to select other 
parameters in proper combination to insure adequate performance. Such 
selection must be compatible with physical realizability at a "reasonable" 
cOSi. 

It has been previously determined that dynamic behavior of a 
ternary system is determined by two non-dimensional quantities: 
sample time, Ve and dead zone, A: Knowledge of the effect of these 
two quantities on system performance, offers the designer a "crutch" 


in design analysis. 


3.2. SELECTION OF NON-DIMENSIONAL SAMPLE TIME 


Non-dimensional sample time is the ratio of sample time to the 
characteristic time of the viscous damped pendulum. It is desirable to 
select system parameters such that limit cycles of higher order are 
prevented. Since the value of non-dimensional sample time has a 
definite effect upon the modes of oscillation possible, it should be 
selected to prevent as many higher order limit cycles as possible. 

Figures C-2 through C-6 show the freedom in position of a simple 
limit cycle about a switching line, plotted as a function of non-dimensional 
sample time. A simple limit cycle is defined to be a limit cycle that is 
complete after two switches between torque modes. Appendix C points 
out that if the freedom in position of a simple limit cycle is negative} 
then the simple limit cycle under consideration cannot exist. In general, 
as the non-dimensional sample time, Uae is increased, more simple 
limit cycles are prevented, It is observed that the 2:2 limit cycle can- 
not be prevented by adjusting the non-dimensional sample time as the 
freedom in position of the limit cycle is positive for all values of sample 
time. Likewise the 1:N limit cycles, where N is any integer, cannot be 
prevented. The 1:N limit cycles are of the lowest mode and are desirable. 

A value of five or greater for non-dimensional sample time will 
prevent all non-desirable simple limit cycles with the single exception 
of the 2:2 limit cycle. It is realized that the simple limit cycles rep- 
resent only a few of the infinity of limit cycles possible. However, one 


would expect that the dynamic behavior of the pendulum would not criange 


Ay = 





greatly for a small change in acceleration, For this reason, it is felt 
that the simple limit cycles are representative of the over-all behavior 


of the pendulum. 


3.3, SELECTION OF NON-DIMENSIONAL DEAD ZONE 


The non-dimensional dead zone must be larger than the non- 
dimensional sample time in order to prevent torquing in a direction to 
aid acceleration. If itis larger than necessary, the bias angle will be 
larger than necessary and errors associated with the bias angle will be 
targer. It would appear that the best value for the non-dimensional 
dead zone is one which is larger than the non-dimensional sample time 


by an infinitesimal amount. 


3.4, SELECTION OF DIMENSIONAL PARAMETERS 


The physical dead zone that is inherent to the error detector is 
a limiting factor in maintaining the null error small. It would appear 
advisable that the dead zone not be made larger artificially, in order 
to implement a ternary pulsed system. 

The key to the proper selection of dimensional parameters lies in 
the dimensional to non-dimensional transformations of Table 2-I. If the 
physical dead zone of the error detector is taken to be the value that will 
transform into a non-dimensional dead zone that is equal to the non- 
dimensional sample time, a system of equations is introduced which 
may be used in the selection of dimensional parameters. The fact that 
the non-dimensional dead zone must be larger than the non-dimensional 
sample time by a small amount, can be taken care of by final adjust- 
ment to the threshold detector. 

Suppose that SUN is the dead zone expected from the error de- 


tector, Then by the transformations of Table 2-1: 


ACS 
ON = (3,1) 
MJ 
and 
Ct. 
T. = ae P (3. 2) 





By desire T. and OA are set equal 


aye (3, 3) 


The applied torque, M, must be sufficient to overcome the reaction 


of the maximum expected acceleration. Therefore: 


Wile a (3, 4) 


It o max? Ans 


fixed by the designer to satisfy the above four equations. The designer 


and te are specified, six parameters remain to be 


may select a value for T. to give satisfactory dynamic performance, 


leaving five parameters to be adjusted to satisfy the four equations. 


3.9, EXAMPLE 


Suppose it is desired to design a ternary pulsed pendulous ac- 
celerometer with a maximum acceleration capability of Loe cm/sec" and 
with a sample time of 500 microseconds. The velocity quantum is the 
product of these two quantities or 5 cm/sec. A value of five is selected 
for non-dimensional sample time so that all undesirable simple limit 
cycles are prevented with the exception of the 2:2 limit cycle. Suppose 
also that the dead zone expected from the error detector is five micro- 


radians. Substituting these values into Equations (3.1) through (3. 4): 


OP aN aee. (3.5) 
M 
Se Woeveee (3,6) 
i 
M. 107 ey a . (aay 
Pp 


If Pis arbitrarily selected to be one gm-cm., Equations (3. 5) 


through (3.7) can then be solved for the remaining unknowns: 


melieye 





M = 107 dyne-cm 


J = 100 gm-cm~ 


and C= 10° dy me aeml— sec, 

Figures 3-1 through 3-ll show the steady state phase plane response 
of such a system for various values of input acceleration. Figures 3-l2 
through 3-15 show the null error versus time response for the phase 
plane responses of Figures 3-1 through 3-ll. 

The phase plane plots show that the pendulum is in its lowest mode 
of oscillation except when the acceleration is in the vicinity of one half 
the maximum. This could be expected in that the 2:2 limit cycle has not 
been eliminated. If the 2:2 limit cycle does not produce excessive error 
in the computer, or if the accelerometer is not expected to operate in the 
vicinity of one half the maximum, then the system is satisfactory in that 


sufficient higher order limit cycles have been eliminated. 


3.6, LIMITATIONS TO THE METHOD 


The authors admit that the previous example has been somewhat 
rigged" so that realistic answers could be obtained for dimensional 
parameters. The specifications could be such that Equations (3. 1) 
through (3. 4) produce dimensional quantities that are not physically 
realizable if a value of five is selected for qe One may have to settle 
for a value much less than five for T.. in which case another method 
will have to be employed to maintain the limit cycles in their lowest 


mode. Such methods are the subject of the next chapter. 
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Null Error in Velocity Counts 
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CHAPTER 4 


COMPENSATION OF TERNARY MODE PULSE TORQUED 
PENDULOUS INTEGRATING ACCELEROMETERS 


oa INTRODUCTION 


Torquing only in the direction to Oppose the acceleration reaction 
torque and torque pulsing ina routine of the lowest order are the en- 
visioned benefits of ternary mode pulse torquing. By routine of the 
lowest order is meant torquing with a routine which keeps the area en- 
closed by the limit cycle as small as possible. When the input accel- 
eration-maximum acceleration ratio r, is less than one-half, this will 
be done if the restoring torque pulse is the shortest time possible, which 
is Z> the non-dimensional sampling period. For values of r greater 
than one-half, the non-torquing period should have the shortest possible 
time duration, again T.. 

Any rational value of r can be represented by a fraction with 


integer numerator and denominator. 





we lS ; (4.1) 
K+ L 
where K is the total number of torquing sample periods per counting 
cycle, and L is the total number non-torquing sample periods per 
counting cycle. A counting cycle is of K + L sample periods duration 
and is tne smallest number of sample periods which accurately rep- 
resents r. Both K and L can be the number of sampling periods in 


more than one torque or non-torque periods. Or 


K = K, + Ky + Kg + SD Oh ot (4, 2) 
where K); Ko, K., -+++ are the number of sampling periods in individuai 
torque periods, and | 

L=L,+blyg+hl2 + oa eee seqenene (4, 3) 
where L,, Lo, L., ---- are the number of sampling periods in individual 
non-torque periods. When r is less than one-half the lowest order counting 
cycle can be produced if K, Ko; K,, -+++ equal one. Likewise, whenr is 


greater than one-half, the lowest order counting cycle can be produced if 
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L,, Lo, Le, eee equal one, 

The above criteria will be the goal in compensating pulse torqued 
pendulous accelerometers. In achieving this goal errors other than 
the error due to higher order counting cycles will be introduced. These 
new errors should be smaller than the errors produced by the larger 
limit or counting cycles, or they must be of a nature such that the 
total velocity stored in the computer can be corrected for the errors 
introduced by the compensation, 

The larger order limit cycles found in sample data bang-bang 
servos are the result of lag introduced in the system by discrete sam- 
pling rather than continuous sampling. Thus, lead must be introduced 
into the error circuit. In nonlinear servo theory, on the error-error 
rate phase-plane, a servo with lead in the error circuit and a dead 
zone has its switching lines represented in the manner shown in Fig- 
ure 4-l. In the case of sample data bang-bang servos, switching from 
one torque state to another occurs at the first sampling instant after 
crossing the switch line. 

In servo mechanisms in which the error magnitude and sense is 
represented in an analog manner by a voltage, lead compensation can 
be inserted by active or passive electronic networks. In pulse torqued 
integrating accelerometers null error is analoged by the voltage output 
of the microsyn. Representative sampling frequencies are one to five 
kilocycles (1000 to 5000 c.p.s.). Maximum frequencies at which micro- 
syns have been found to be excited are also one to five kilocycles. Phase 
distortion increases as the exciting frequency is increased. The error 
analog is of poor quality. In existing binary pulse torqued pendulous 
accelerometers, the microsyn null position error detector is used 
only to determine the sense of the error. Hence, because of poor error 
analog quality, compensation of the error analog signal is not practical. 

What is available in the system is the error sense sampled at 
equal periods. This suggests compensation using logic techniques. 

To simulate the lead switch lines shown in Figure 4-l, error sense, 
magnitude, and rate must be obtainable. The best magnitude resolution 
which can be expected from a microsyn type position error detector is 
about the width of the microsyn dead zone. Error rate is not readily 
available. The most which can be obtained readily from microsyn null 


error detectors at the desired sampling frequency is pendulum position 
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Figure 4-1. Phase Plane Representation of a Second Order 
Position Servo With Dead Zone and Lead. 
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iri zones, These zones are the dead zone and magnitude zones of 
detectable sense, 

The above zones can be determined and controlled in width with 
limited ability by passing the microsyn output signal through a phase 
detector to determine sense, Then the signal can be applied to biased 
level detectors to set zone magnitude, 

Figure 4-2 shows the phase plane divided into five zones. 
Following are logic techniques in using these zones to insert lead 


compensation into the system. 


4,2, ALTERNATING PULSE ZONE COMPENSATION 

Lead compensation can be interpreted as a means of anticipating 
required actions and carrying them out sooner than a non-compensated 
system would. Figure 4-1 shows that the magnitude of the lead effect 
should be proportional to the magnitude of the velocity. This is not pos- 
sible with available error detectors at the information rate required. 
Where our system does not allow a pure lead effect, the lead producing 
zone can produce a quasi lead effect. 

A zone of sampling periods of alternating torquing and not torquing 
periods will anticipate both the change from the non-torquing dead zone 
to the continuous torquing zone and vice versa. Considering positive 
input accelerations, it is obvious for values of input acceleration - 
maximum acceleration ratio, r, greater than 0.5 that more torquing 
periods are required in a counting cycle than non-torquing periods. 
Hence, @ will increase until switching is about the boundary between 
the continuous and alternating zones, the-@ line, Similarly for values 
of r less than 0.5 the switching will be about the boundary between the 
alternating zone and the dead zone, the -¢pline. 

Previously the non-dimensional equations of motion of the pendulum 


were determined to be 


=(6-rrte -l+(-e )b +6, (4. 4) 


T 


Ae ee (4.5) 


=6.9- 





"SsouOZ SUTONpOIg pea] oy} SuTMOYS SsUuoZ OU] aueTd asseud at} Jo UOTSTAIG 


sUu0Z7 
anbIo, 9U07 
SATICSIN suTONpoOId 

snonutjuoz peayT 


: 
| 
| 
| 
| 
| 


| 
| 
| 


IU0ZF 


peod 


| aUu07 

sutonpold 
| peaT 
| 


"C-h einsi A 





IU07 
anbso0 
9ATITS Of 

snonutjuog 


~70- 





Except for the situation in which a very quick change in the accelera- 
tion has reversed the sense of the acceleration such that the accelera- 
tion reaction torque and the applied error reduction torque are in the 
Same direction, the magnitude of maximum ¢ in the alternating torquing 
zone is seen from equation (4.2) to be (6 - r). One-half is the maxi- 
mum value of r when switching from the dD line. Since when moving 
from the dead zone into the alternating torque zone 6 equals zero, the 
maximum b associated with switching about the “on line in the alter- 
nating torque zone is 0.5. Similarly, since 0.5 is the minimum value 
of r and 6 = 1 when switching about the “Or, line and moving from the 
continuous torque zone into the alternating torque zone, the maximum 
magnitude of b associated with switching about the “Or, line in the 
alternating torque zone is 0.5. 

If b is greater than 0.5 in the alternating torque zone the 
switching locus is moving from one boundary to the other because of 
a Change in r from larger than 0.5 to smaller than 0.5, or vice versa. 
In this case switching in the alternating zone is not necessary. If $ is 
less than 0.5 switching is mandatory. This will be assured if the 
width of the alternating zone is made at least equal to the maximum b 
desired to be "caught" in the alternating zone times the non-dimensional 
sampling period T.. The minimum value of Or, = oD acceptable is 
0.5 T.. 

The width of the dead zone must be adjusted so that, for very 
small values of r, the pendulum does not pass into the alternating zone 
where the magnitude of (6 - r) is greater thanl. Refer to the conditions 
shown in Figure 4-3 for the limiting case of small r, that is, r approxi- 
mately equal to zero, %, equals “Op: b, approximately equals zero. 
~o is less than top: %o is the point of furthest penetration into the 
positive @ plane before the pendulum is driven back to the negative ¢ 
plane. By equations (4. 4) and (4.5), with a torque pulse of magnitude 
(6 - r) approximately equal 1, the position one sampling period later is 


= 10 


me ee (4, 6) 


oD 


$= 1-e ° (4, 7) 
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Figure 4-3, Phease Plane Condition for Determining Minimum Width 
of the Dead Zone in a Compensated Ternary Pulsed In- 
tegrating Pendulous Accelerometer. 
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When the pendulum is coasting, that is, 6 equals zero and r equals 
zero, and the phase plane path has not crossed the zero velocity axis, 


the motion can be described by 
Oye) (4, 8) 


Hence, when @ equals zero at 5 


bo =O + (4, 9) 
-T = 

do=T te -+-Il-¢,t+l-e ~ (4. 10) 

do = T. - $p: (4, 11) 


To keep @ in the dead zone, 
(+d) - (-$p) > by - (dp), (4, 12) 
(+67) - (-@p) > T., (4, 13) 


or, the non-dimensional magnitude of the minimum width of the dead 
zone is equal to the magnitude of the non-dimensional sampling 
period, T.. 

An elemental switching logic is explained and designed in 
section D. 2 of Appendix D. A symbolic design is shown in Figure D-4, 

Figures (4-4), (4-5), and (4-6) show the steady state counting 
cycles for the alternating pulse zone compensation for values of r in 
one-tenth increments from 0,1 to 0.9. These values were chosen 
because the counting cycles are relatively simple yet they exhibit the 
problems found in the technique. Where practical, two counting 
cycles are shown and they indicate the boundaries of the area in the 
phase plane in which the limit cycles can exist. 

For the cases of r equal 0.3 and 0.7 only one limit cycle is 
shown. Inspection of the cycles show they have a very limited freedom 
of position, Pr The trajectory of the error on the phase plane will 


adjust itself and hunt for the limit cycle. In Figure (4-4) the area in 
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Figure 4-6, Steady State Limit Cycles for Alternating Pulse Zone 
Compensat:on Me =1,0, 
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wich the trajectory for r equal 0.3 hunts, after large transient effects 
have died out, is shown with dashed lines. One can see that the over- 
all counting error during the hunting period is about that of the steady 
state limit cycle. 

For values of r less than 0.5 the counting cycles are found to be 
in the lowesi order. For values of r greater than 0.5 requiring a two 
period tordue pulse the counting cycles are not in their lowest order. 

r equal 0.6 requires a torque 2 - not torque l, torque 3 - not torque 1 
counting cycle for lowest order instead of 3-1, 1-1, 1-1, 1-1 counting 
cycle. r equal 0.7 requires 3-1, 2-1, 2-1 counting cycle instead of 
o-l, 1-1, 1-1 counting cycle. 

In order to have a two period pulse in the alternating pulse zone 
compensation, the pulse must start in the continuous positive torque 
zone and move into lead compensation zone before the end of one sam- 
pling period in time. For different values of r, meeting the above con- 
dition depends on the value of TT. the non-dimensional sampling time, 
or the ratio of the sampling period to the time constant of the pendulum. 

Figures 4-7 through 4-12 show the steady state counting cycles 
for alternating pulse zone compensation with T equal 5.0. Those 
values of r which were in their lowest order at T. equal l.0O havea 
larger freedom of position at a equal 5.0. For r equal 0.6, 

Figure 4-10, the lowest order counting cycle can now be assumed with 
some freedom of position, Pr For r equal 0.7, Figure 4-9, the lowest 
rder counting cycle is possible, but the system is required to switch 
at almost a point on the phase plane. The system will "hunt" for this 

point in the 3-1, 1-1, 3-1 order after large transient effects have 
died out. 

Figure 4-11, for r equal 0.88 and tee equal 5.0, shows the 
counting cycle is in its lowestorder; 8-1, 7-1, 7-1. At ue equal 1.0 
the Pe is approximately one-fifth that shown in Figure 4-12, again re- 
quiring more "hunting" to get in a lowest order counting cycle. For 
larger values of r which do not require a two sample period pulse in 
their lowest counting order, the lowest order is possible for small 
values of T., but increasing 7 reduces the amount of "hunting" re- 


quired to assume the lowest order counting cycle. 
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Figure 4-7, Steady State Limit Cycles for Alternating Puise Zone 
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Figure 4-8. Steady State Limit Cycle for Alternating Pulse Zone 
Compensattor T. = 9, OF 
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Figure 4-9. Steady State Limit Cycles for Alternating Pulse Zone 
Compensation T_ = 95.9. 
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Figure 4-li. Steady State Limit Cycles for Alternating Pulse Zone 
Compensation Jee = 9,0. 
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4,3. INVERSE ORDER ALTERNATING PULSE ZONE COMPENSATION 


in the alternating pulse zone compensation enough lead effect 
was introduced as the position entered the lead compensation zone from 
the dead zone. However, on entering the lead compensation zone from 
the continuous torque zone the system continues torquing one sample 
pericd after the first sample in the lead compensation zone before it 
switched to non-torquing. This removed most of the lead effect. By 
changing the switching logic so that when a sample is in the lead com- 
pensation zone the action during the period following the sampling 
instant is opposite the action carried out the previous period, that is 
non-torquing if torqued previous period and vice-versa, an inverse 
order alternating pulse zone is set up. This will produce the desired 
lead effect in both the increasing and decreasing @ direction. 

For the same reasons as in the alternating pulse zone compen- 
sation, the width of the dead zone must be greater than T and the width 
of the lead compensation zone greater than 0.5 ie 

The switching logic for the inverse order alternating pulse zone 
compensation is developed in section D-3 of Appendix D and shown in 
the diagram of Figure D-5. 

The lead effect when switching about the oD lines was not 
changed, The phase plane results for values of r from 0 to 0.5 are 
the same as for the alternating pulse zone compensation shown in 
Figures 4-4 and 4-5. The phase plane results for Te equal l andr 
greater than 0.5 are shown in Figure 4-l2 and 4-13, The resulting 
lowest order counting cycles are shown. Also Figure 4-13 shows the 
expected result, that if Tes is increased the counting cycles remain 
in their lowest order. 

If the compensation problem is viewed from the aspect of not 
allowing higher order limit cycles by parameter variation as de- 
veloped in Chapter 3, the same results are obtained as in the inverse 
order alternating pulse zone compensation except as to the width of 
the compensation zone. 

Referring to Figure C-l, if the portion of the switching line for 


ositive angular rate is moved to the left of its present position 
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Figure 4-12, Steady State Limit Cycles for Inverse Order 
Alternating Pulse Zone Compensation iT. oe el 
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Figure 4-13. Steady Stute Limit Cycles for Inverse Order 
Alternating Pulse Zone Compensation. 
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the N= K - 1 sample will be to the right of the switch line and the 
system wil! stop torquing. Similarly if the portion of the switching 
line for negative angular rate is moved to the right of its present 
position, the N = K + L - 1 sample will be to the left of the switch 
line and the system will start to torque. The action described has 
the effect of making the freedom in positon, Pes negative for higher 
order limit cycles. 

For the iowest order limit cycles, and r greater than 0.5, the 
system may non-torque only one sample period, T., ata time. The 
positive @ portion of the switching line can best be instrumented by 
the positive portion of the Or, Switching line of the inverse order 
alternating compensation zone. By similar reasoning, the negative 
$ portion of the switching line can best be instrumented by negative 
portion of the oD switching line of the inverse order alternating com- 
pensation zone. Subsequent samples still in the compensation zone 
will cause alternating torquing. 

From Figure C-2 we see if the switching line separation is made 
greater than 0. 363 for ie equal 1.0 the 2:2 limit cycle is eliminated. 
Figures C-3 through C-6 show if Pe eliminates the 2:2 limit cycle, 
all higher limit cycles are eliminated. ly, “ on! effects g,. Inspec- 
tion of Figure C-2 shows that as T increases the required P does 
not increase linearly. In fact it does not exceed 1.0. 

The logic required is the logic for the inverse order alter- 
nating torque zone compensation. Viewing the compensation problem 
from the parameter variation viewpoint resulted in the same com- 
pensation logic but with the minimum width of the dead zone decreased 
from 0.95 ie to 0. 363 oe 

The result of decreasing the width of the lead compensation zone 
is to decrease the error bias and cross coupling effects. But for very 
large transient situations the time required to reach a steady state 
limit cycle will be longer. 

The phase plane trajectories for T. equal one are similar to 
Figures 4-12 and 4-13 except that the “Or, line is at -0. 363 instead of 


-0.5 and thus the error bias is reduced by 0.137. 
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4.4, LOGICALLY FORCED LOWEST ORDER COMPENSATION. 


In this compensation, switching logic is devised to allow the sys- 
tem to torque only one period at a time if r is less than 0.9 and to not 
torque only one period at a time if ris greater than 0.5. Also limit 
cycle switching will be limited to be about the %D line in order to re- 
duce cross coupling errors. If ris greater than 0.5 and the system 
only torques one period at a time the limit cycle position will move 
toward greater error. If ris less than 0.5 and the system will not 
torque only one period ata time, the limit cycle position will move 
fo crmaliler error. 

Based on the described movement, when it is desired to switch 
mrom logic for r less than 0.5 to logic fer r sreater than 0 osu 
position will be moving toward greater error. If a zone or boundary 
line called tor or Or, aS appropriate, is set up beyond the greatest 
error of counting cycles for r less than 0.5, whenr is greater than 
0.5 the position will increase until it passes this line. When a sam- 
ple with @ greater than Or is obtained the logic can be changed to 
torque as required when out cf the dead zone, but not torque less 
than one sampling period, Ts: ata time. The system will then torque 
until the position returns to the dead zone and then switching will be 
about the D line. This logic system requires the position to return 
out of the dead zone in only one period of not torquing. To insure this 
for all combinations of dead zone penetration and all values of r 
preater than 0.5 it is found that a must be at least 5. 

The dead zone width requirements are the same as before, 
(+95) - (-¢) greater than T.- 

Then with as equal 5 or greater, r will be less than 0,5 when 
two consecutive samples are found to be in the dead zone. This can 
be used to change the switching logic when r becomes less than 0. 5. 
For consecutive samples in the zone between Pr, and %p» the system 
will torque alternating periods. 

The logic for this compensation is developed in section D-4 and 
shown in Figure D-6. 

Figures 4-14 through 4-18 show the counting cycle for the 


logically forced lowest order compensation to be in the lowest order. 
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Figure 4-14. Steady State Limit Cycles for Logically Forced 
Lowest Order, ee = 3, 0) 
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Figure 4-16. Steady State Limit Cycles for Logically Forced 
Lowest Order, T. = OG 
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Figure 4-17. Steady State Limit Cycles for Logically Forced 
Lowest Order, di =u9. U8 
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Figure 4-18, Steady State Limit Cycles for Logically Forced 
Lowest Order, Te = 9.0. 
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In fact for values of r other than in the 0.5 vicinity the counting cycles 
are exactly the same as in the non-compensated system with T, equal 


09.90. Only in the vicinity of 0.5 is the counting cycle improved. 


4.9. COMPENSATION RESULTS 


On each of the figures showing limit cycles the zero error line 
is shown, This allows one ata glance to get a relative estimate of 
tne error size. In Figure 4-19 through 4-24 the values of the error 
during the limit cycles are shown. The errors were computed using 
the previously developed 

E=o+¢. (4, 14) 


E is the non-dimensional error expressed in counts of velocity quanta. 
The largest error allowed by the freedom of position is shown. The 
following conclusion can be drawn from the error plots and the phase 
plane trajectories. 

In the alternating pulse zone compensation T. must be in- 
creased to 5 to insure the counting cycles are in their lowest order. 
The only performance improvement is in the vicinity of r equal 0.5. 
The hardware restrictions placed on the system by the uncompensated 
system due to a large se has not been removed. In addition, for values 
of r greater than 0.5 the error bias has doubled. The more compli- 
cated error detection and switching logic is a relatively high price to 
pay for doubtful performance improvement. 

The logically forced lowest order compensation also requires 
T. greater than 5. It has the same error bias as the uncompensated 
system. Its counting cycles are the same as those of the uncom- 
pensated system with i equal 5 except when r-is in the vicinity of 
0.5. The logically forced lowest order compensation does produce 
the lowest order counting cycles for all values of r. 

The inverse order alternating pulse zone compensation pro- 
duces the lowest order counting cycle for all values of r with T, equal 
one or greater. For r equal 0.5 and greater the error bias is greater 
than in an uncompensated system due to the switching about the or 
line. The switching logic shown in Appendix D for the inverse 


alternating pulse zone compensation is relatively simple. 
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Figure 4-19, Error Cycles for Compensated Systems With Te = 1,0 
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Figure 4-21. Error Cycles for Compensated Systems With T. = ie 
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Figure 4-22, Error in Compensated Systems With t. = 9,0 
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The velocity quanta size, the maximum acceleration to be meas- 
ured, and the computer clock frequency usually are fixed by the appli- 
cation of the pulsed pendulous accelerometer. The dead zone of the 
error detector is fixed in minimum size. These fixed parameters, 
with the additional requirement that the width of the dead zone be 
greater than one T., place stringent demands on the ratio of the pen- 
dulum moment of inertia J to the damping constant C. Smaller mini- 
mum required T eases the restriction on the J to C ratio. 

With the non-dimensional width of the dead zone equal to one T., 
the error bias due to switching about the DH line is one-half a velocity 
quanta and that due to switching about the Or line is one velocity 
quanta. Sections D.5 and D.6 of Appendix D develop logic to correct 
the velocity count for the error bias due to switching about particular 
switching lines. This correction makes it practical to make the dead 
zone more than one ite wide in non-dimensional magnitude. This will 
ease the restrictions put on the J to C ratio. Correcting the velocity 
count for the error bias does not decrease the physical angle the pen- 
dulum makes with its null position, It will not decrease cross coupling 
effects. 

With proper error bias correction and using the inverse order 
alternating pulse zone compensation, the counting error can be kept 
within one quanta. 

In describing the limitations of the error detector the conserva- 
tive statement was made that the best magnitude resolution to be ex- 
pected was about equal to the width of the dead zone. This assumption 
was based on the high noise to signal ratio due to a relatively low error 
information rate. If quadrature effects about the microsyn null posi- 
tion are the major factor in creating the error detector's dead zone 
size, better resolution away from the dead zone can be expected. 

Better resolution away from the dead zone will allow the non- 
dimensional dead zone to be several T. wide, but the distance be- 
tween ¢, and oD still only needs to be 0.5 T.. $7, can then be moved 
in towards the dead zone decreasing cross coupling effects. This will 
require error count correction because the dead zone will become 


the several counts wide. Also, because the dead zone is several 


-100- 





counts wide the random initial error increases. This could be pre- 
vented by electrically restraining the pendulum until the count is 
initiated. This would require a special counting error bias initial 


correction 


4.6. COMPENSATION SUMMARY 


Three logic compensation schemes were devised. The alter- 
nating pulse zone compensation with a small amount of switching logic 
traded larger error bias for keeping the counting cycle in its lowest 
mode for all values of r. It required T. equal to five or more placing 
design restrictions on the hardware. The logically forced lowest 
order did not increase the error bias in making all counting cycles 
the lowest order, but it again requires iy equal five or more and uses 
more complicated switching logic, The inverse order alternating pulse 
zone cOmpensation puts the counting cycles in their lowest order and 
only requires es equal to or greater thanone. The switching logic is 
Simple. For values of r greater than 0.5 the error bias is increased, 
Velocity count correction schemes were devised to correct the counted 
velocity for the existing error bias. The inverse order alternating 
pulse zone compensation is the most practical because it does not 
place stringent moment of inertia - damping coefficient ratio re- 


strictions on the pendulum unit. 
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CHAPTER 95 


CONCLUSIONS AND RECOMMENDATIONS 


9.1. CONCLUSIONS 


The theoretical feasibility of ternary pulse torquing a pendulous 
accelerometer has been shown, The inherent greater stability of a 
ternary system as compared to a binary system results in smaller 
limit cycles associated with smaller velocity quanta counting error, 

By adjusting the ternary pulsed accelerometer system's com- 
ponents so as to make the non-dimensional sample time, T 5 5,0 or 
greater, the lowest order limit cycles can be assured. When practical 
considerations do not allow adequate adjustment, switching logic com- 
pensation can be used to insure the lowest order limit cycles with Ts 
equal or greater than1.0. In compensating the ternary mode pulsed 
pendulous accelerometer increased error bias is traded for smaller 
limit cycles, If the error bias is too large, the velocity quanta count 
in the computer can be corrected for the error bias. Another result 
of the compensation is increased cross-coupling effects. These can- 
not be as easily corrected in the velocity count and thus must remain 
small, 

Ternary mode of pulse torquing pendulous integrating accel- 
erometers allows increased accelerometer performance to be "bought" 


with increased system complexity. 


9.2. COMBINATION SYSTEM 


Experience gained in studying the ternary mode pulse torqued 
pendulous integrating accelerometer leads to the consideration of a 
combination system. Ternary in mode of error detection, that is, the 
system senses clockwise error, dead’ zone, and counter-clockwise 
error; and binary in mode of torquing, that is, torque clockwise or 
counter-clockwise. The ternary mode in error detection would be 
used to effect lead in the torque switching. The sense of torquing would 
be changed on entering the dead zone, Investigation to determine if 
enough lead would be introduced to insure lowest order limit cycles is 


recommended, 
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9.3, RECOMMENDATIONS FOR FURTHER INVESTIGATION 


The hardware problems specific to ternary pulse torquing have 
not been investigated, Investigation to determine the feasibility of 
having a system produce a reproducible torque pulse in two senses and 
also not produce it for specific periods needs to be carried out. An 
energy sink might be the solution except in the case when the ternary 
mode is used to measure small accelerations at little expenditure of 
energy. The environmental effect of a varying rate of energy expen- 
diture by the torquer should be studied. The exact performance of 
error detectors at high sampling rates both in the dead zone region and 
away from the dead zone is needed to insure the feasibility of the com- 
pensation schemes, 

Investigation in the use of ternary mode pulse torquing gyros, 


gimbals, and other inertial instruments is recommended. 
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APPENDIX A 


NOTATION 


Dead zone of error detector 

Angular Position of the pendulum with respect to the null 
position 

Component of acceleration along the input axis of the 
accelerometer 

Component of acceleration along the Reference Axis (or null 
axis) of the accelerometer 

Coefficient of viscous friction 

Moment of inertia of the pendulum 

Time 

Time between samples 

Torque Generator torque 

Magnitude of Torque Generator torque 

Pendulosity of the Pendulum 

Velocity increase 

Indicated velocity 

Total error in velocity indication 

Digital error 

Null error 

Cross-coupling error 

Torquing mode, +l, 0, or -l. 

Non-dimensional null angle 

Non-dimensional time 

dg/dr 

Non-dimensional dead zone 

Non-dimensional acceleration along input axis 
Non-dimensional acceleration along reference axis 
Non-dimensional velocity 

Non-dimensional velocity error 

Norn-dimensional null error 


Non-dimensional cross-coupling error 
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Non-dimensional digital error 


E 
d 
ie Non-dimensional sampling period 
Pr Freedom in position of a simple limit cycle about a switching line 
> Equal by definition 
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APPENDIX B 


DERIVATION OF ERROR DETECTOR DEAD ZONE REQUIRED TO 
PREVENT A LIMIT CYCLE WHEN INPUT ACCELERATION IS ZERO 


Consider a limit cycle of K time intervals in the positive torque 
mode followed by L time intervals in the zero torque mode, K time in- 
tervals in the negative torque mode and L time intervals in the zero 
torque mode, The limit cycle-must be symmetric in that the sum of the 
torque pulses over a limit cycle must be zero. Such a limit cycle is 
shown in Figure B-l. Two requirements must be met in order that the 
limit cycle can exist: 


eK+L-1) ~ 9(K-1) = 9A a 


Tieuuey Piel) A (B-2) 


pT 


Now suppose that switching line 1 is held fixed and switching line 2 
is moved to the left until Condition (B-1) cannot be met. Switching will 
occur prematurely and K will be reduced by one. Since the distance 
traveled in the torque mode is reduced, the distance traveled in the no 
torque mode must increase so that Condition (B-l) can be met for the 
new value of K. Thus L must increase. 

Switching line 2 can be moved to the left until K is reduced to one 
and L is increased to infinity. As L approaches infinity the velocity at 
N = K+L must approach zero since the pendulum has been damped 
wihout being forced for an infinite amount of time. The width of the 
limit cycle under these conditions can be determined from Equations 
(2-7) and (2-8): 

-T =f 
bea)“ to* [Te --1 +¢1-e * (B-3) 
or 


9(«) = %o 7 ae (B-4) 


Now, if switching line 2 is mOVed to the left such that PA is greater 


than TS: then K must be reduced to zero and the limit cycle cannot exist. 
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Figure B-l. Limit Cycle With Zero Acceleration. 
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Thus the condition necessary to prevent a limit cycle is: 


On aie (B-5) 
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APE ENDESS 


DERIVATION OF THE FREEDOM IN POSITION OF A SIMPLE 
LIMIT CYCLE ABOUT THE SWITCHING LINE, SHOWING THE 
CONDITION NECESSARY FOR A SIMPLE LIMIT CYCLE TO EXIST 


A simple limit cycle is defined as a limit cycle which is com- 
plete after two switches. For the purposes of this derivation it is 
assumed that the zero switch zone is of sufficient width to insure that 
the pendulum in steady state remains in the positive and zero switch 
zones for positive acceleration and in the negative and zero switch 
zones for negative acceleration. A typical limit cycle consisting of 
K sample intervals in a torque mode followed by L sample intervals 
in the zero mode is shown in Figure C-l. 


The freedom in position of the limit cycle is defined as follows: 


A 
JEST beh) © = Th (c.0) 


It is evident from Figure C-1l that the complete limit cycle may 
be shifted to the right or left over a total range of OF without disturbing 
the switching sequence. It is also evident that if be is negative, the 
switching sequence must be altered and the limit cycle described 
cannot exist. 

The velocity summing capability of the system depends upon the 
relation 

K 


(C2) 
|e en 





r = 


Equation (C. 2) is without error since the change in null error Ey 
over a limit cycle is zero, Since K and L are integers, the input ac- 
celeration must be the ratio of two integers if the limit cycle is to exist. 


The equation of motion of the pendulum is: 


$+¢6=(6 -7r) (C. 3) 


where 6 is the torquing mode. Equation (C. 3) is linear between torque 
switches and may be integrated as such with results being valid be- 


tween switches, 
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Figure C-l. Simple Limit Cycle, 
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, mn - 7) -(T-T.) 
_ AE se ee oie | eae : (C. 4) 
: Be) -(T-T.) 
$=, +@ \l-e + (0"= ) Gaeta — (C. 5) 
i i 


A solution for O(KEL) in terms of 6, is obtained using Equation (C. 4). 


$ -(K+L)T -L,T -(K+L)T -(K+L)T 
fen) %o © e+e e - e mene eer (C.6) 


where Te is the non-dimensional sampling interval. 


Since a limit cycle is complete after K plus L sample intervals, 





%o 3 %K+L) (C, 7) 
Substituting (C. 7) into (C.6) and solving for bo: 
-LT 
; e eal 
meee (Comn:)) 
(K+L)T . 
l-e 
Substituting (C. 2) into (C. 8) 
-LT. 
ee (c. 9) 
AK+L)T. K+L 
l-e 


By using (C. 4), (C. 5), and (C. 9) an appropriate number of times 





ers. Ricleaiiek. - (K-DT, LT(K-1) 
ee ee Oe 
SU Selly) 20 K+L 
l-e 
and 
$ =o + is ac: iu 28 Bile @2) Ks on 
l-e 
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Substituting (C.10) and (C.11) into (C. 1) 


(K+L-1)T. sce (L-1)T. 


¢.==8 et tae Ue ee (Cul) 
AK+L)T. K+L 





l-e 


Figures C-2 through C-6 show De as a function of the non-dimensional 


Sampling period for various values of K and L. 
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Freedom in Position of Simple Limit Cycles About a 
Switching Line; 1:N Limit Cycles. 
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Figure C-3., 


VAS De 





Freedom in Position of Simple Limit Cycles About a 
Switching Line; 2:N Limit Cycles. 
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Figure C-4. 





Freedom in Position of Simple Limit Cycles About a 
Switching Line; 3:N Limit Cycles. 
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Figure C-5. Freedom in Position of Simple Limit Cycles About a 
Switching Line; 4:N Limit Cycles. 





Figure C-6. Freedom in Position of Simple Limit Cycles About a 
Switching Line; 5:N Limit Cycles. 
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APPENDIX D 


SWITCHING LOGIC DESIGN 


D.l. INTRODUCTION AND DEFINITIONS 

The object of this appendix is to present the relative complexity 
of the different switching logics required. 

The logic design will be at the logical element level, A logical 
element is a computer element which performs a nontrivial logic 
function; that is, some input or combination of inputs produces 
something new as an output. 

In this presentation of switching logic design elementary know- 
ledge of switching design is presumed, Chapter 17 of Reference 6 is 
an elementary presentation of logical design and was used in de- 
veloping the switching logic presented here. Switching functions will 
be in Boolean algebra form. 

Figure D-l defines symbols and characteristics of the logical 
elements used, The delay repeat is not a standard element. Figure 
D-2 is an elemental layout of the delay repeat. 

Figure D-3 shows the definition of the error level signals. The 
error detector's output will be zero in the dead zone; but D, the dead 
zone position signal, can be computed, since it is the absence of the 
other error level signals. 

Table D-1 defines the letter symbols applicable to this appendix 


only used to define quantities used in the logic design. 


D.2. SWITCHING LOGIC FOR ALTERNATING PULSE ZONE COMPENSATION 
The logic for switching in the alternating pulse zone compensation 
can be represented by Boolean algebra as follows: 
Always torque when in the torque override zones, 
P=F (D. 1) 
N=J (D, 2) 


When in the lead compensation zone, torque One T do not torque 


next T. if still in the lead compensation zone, and then repeat the cycle 
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Figure D-l. Symbols and Characteristics of Logical Elements. 
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Output 


Figure D-2, Elemental Design of the 
Delay Repeat. 
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W KK x x < 


TABLED -!1 
DEFINITION OF LETTER SYMBOLS FOR SIGNALS 
USED IN SWITCH LOGIC DERIVATIONS ONLY 


- Signal A at present time instant 
- NotA 
- Value of A at first preceding interval 


- Value A has maintained constant for x preceding periods. 


x may take any integer value, 


- Positive torque 

- Negative torque 

= Error ain dead @zone 

- Error in positive continuous pulse override zone 
- Error in negative continuous pulse override zone 
- Error in positive lead compensation zone 

- Error in negative lead compensation zone 

- Positive torque pulse as a result of G 

- Negative torque pulse as a result of H 

- Intermediate signal in computing K 

- Intermediate signal in computing L 


- Mean error bias correction magnitude when switching about 


@ equal (- or +) oD 


- Mean error bias correction magnitude less U/2 when switching 


about @ equal (- or +) Or 


- Negative error sense 

- Permit positive torquing for more than one T. 
- Positive torque for not more than one i 

- Permit negative torquing for more than one it 
- Negative torque for not more than one ae 


- Error has been in the dead zone this sampling instant 


and the one previous 
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as long as in the lead compensation zone, 














K=G-K (es) 

ee (D. 4) 
Combining 

P=-F+K (D. 4) 

P=-Fric- kit! (D. 5) 

N =e) zee (D. 6) 

Neon Ee ee (D. 7) 


The switching logic described by equations (D. 5) and (D.7) is 


shown in Figure D-4, 


D.3. SWITCHING LOGIC DESIGN FOR INVERSE ORDER ALTER- 
NATING PULSE ZONE COMPENSATION 
The switching logic for the inverse order alternating pulse zone 
compensation is like the switching logic for the alternating pulse zone 
compensation except that when in the lead compensation zone the 
opposite action of that taken the preceding period, no matter which 
zone determined the action, is desired. This is done by changing 


equation (D. 3) to 


K=G: P (D. 8) 








and equation (D, 4) to 
ptt 
Lg? Hee Nee (D. 9) 
Then 
Perce poe (D.10) , 
Tt+l 


N=J+tH: N (D, 11) 


The switching logic described by equations (D.10) and (D. 11) 


is shown in Figure D-5, 
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Figure D-4, Switching Logic For Alternating Pulse Zone Compensation. 
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Figure D-5. Switching Logic for Inverse Order Alternating 
Pulse Zone Compensation. 
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D.4. LOGICALLY FORCED LOWEST ORDER COMPENSATION 
SWITCHING LOGIC 

In the logically forced lowest order compensation, for values 
of r less than 0.5 the system is allowed to torque for only one sam- 
ple period, livag at a time,and for values of r greater than 0.5 the 
system is allowed to not torque only one sample period, To, ata 
Hime. The lopic will be chanred trom, multiple we positive torquing, 
X, to single T. positive torquing, X; or from multiple i negative 
Porauing, Y, tO single - negative torquing, Y; when the error 
Signal remains in the dead zone for two consecutive sample periods. 


This is represented by Booloean algebra by 


T+ 


B=D-.- D ; (a2) 
reset x 

xB. (D. 13) 
and reset Y 

ee ae (D. 14) 


The logic will be changed from single = torquing, X or Y, 
to multiple T, torquing, X or Y respectively, when the error signal 


passes the +O or “Oy, line respectively in an increasing error 


direction, 
Set X 
5 
x ies (D. 15) 
Set Y 
& 
We yo (D. 16) 


The torque switching can be summarized as follows: 


T+ 


JE) ey MEG) ge 0s (Gl al 


(D. 17) 


io i eee yet (D. 18) 
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The logic described by equations (D,12) through (D.18) for the 


logically forced lowest mode compensation is shown in Figure D-6. 


D.5. ERROR BIAS CORRECTION SWITCHING LOGIC 


Proceeding from positive error to negative error, the change 
in error bias correction is +V when changing from +o, to top,» 
+U when changing from +O to “Op: and +V when changing from 
“9p to “oy: The sign of the change in error bias correction is 
negative proceeding from negative to positive. The magnitude of 
U and V may be adjusted to make the best correction for the mean 
error due to switching about a line. 

When it is necessary to correct for the error only when switching 


about the Or lines, the Boolean algebra equations for the switching are 


al yas 2 ees = ye (D. 19) 
2 

aie ep nse 1, . p? (D. 20) 
2 


Only one quantity may be stored ata time. Figure D-7 shows the 
required switching logic. 

If it is necessary to correct only for the error bias in 
switching about the oD lines as in the uncompensated case or the 
logically forced lowest mode case again a special starting condition 
need not be used unless the pendulum was restrained until counting 
started. This is because the pendulum will be switching about the 
initial switch line when counting starts. +U must be switched when 
entering the F or G zones when -U was last switched. -U must be 


switched when entering the H or L zones when +U was last switched. 


pa (ie 1h) oy aes (D. 21) 


<= si 4 1b) ena es (D. 22) 
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Figure D-6., 





Torque Switching Logic for Logically Forced 
Lowest Order. 
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Figure D-7. 





@,, Switching Line Error Bias Correction 
Switching Logic. 
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-U will reset (+) Pt and +U will reset (-u) PTX, 


Figure D-8 shows the two state oD Switching line error bias 


correction switching logic. 
A four state error bias correction to correct for the error 


bias due to switching about each of four possible switching lines can 


be constructed using a combination of the two preceding error bias 


corrections. 
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Figure D-8, 


+U 
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@p Switching Line Error Bias Correction 
Switching Logic. 
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